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Abstract 

The two-dimensional ferromagnetic anisotropic Ashkin- Teller model is in- 
vestigated through a real-space renormalization-group approach. The crit- 
ical frontier, separating five distinct phases, recover all the known exacts 
results for the square lattice. The correlation length (z^j-) and crossover 
{(p) critical exponents are also calculated. With the only exception of the 
four-state Potts critical point, the entire phase diagram belongs to the Ising 
universality class. 
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I. INTRODUCTION 



The Ashkin- Teller (AT) model was introduced to study cooperative phenomena of 
quaternary alloys on a lattice. Lately it was shown that it could be described in 
a hamiltonian form appropriated for magnetic systems In this representation the 
AT model can be considered as two superposed Ising models described by classical spins 
variables a and r. In addition, the Ising systems are coupled by a four-spin interaction 
term. The isotropic AT (lAT) model corresponds to the case where the two Ising systems 
are identical to each other. 

The phase diagrams for the lAT model are relatively well known in both two and three 
dimensions. Many exact results were obtained for the lAT on a square lattice 0,^]. The 
corresponding phase diagram was determined by real-space renormalization group analysis 
0, mean- field renormalization-group approach Monte-Carlo renormalization group 

and Monte Carlo simulations [rU|,rT|. These results may be compared to the phase 



diagram of the selenium adsorbed on Ni surface, which is a good physical realization for 



the two-dimensional lAT [IT^jITI]- The oxigen ordering in Y Ba2CuzOz may also be un- 
derstood in analogy with the two-dimensional lAT model |T^-|T6| . The three-dimensional 
lAT model was extensively investigated by a number of techniques , showing an enor- 



mous richness of critical behaviour. Mean-field calculations suggest the existence of many 



multicritical points, some of them confirmed by Monte Carlo simulations [|T^. A Cayley 
tree formulation suggests that, for high enough coordination, there may be also an addi- 
tional phase in the antiferromagnetic lAT where a symmetry breaking between the spins 
a and r takes place |T8[ . 



The physical description of the anisotropic AT model (AAT), however, is not so well 



known. Using exact duality relations Wu and Lin [|T9| determined the general topology of 
the phase diagram for the AAT on a square lattice. Their results were lately supported 
by a Migdal-Kadanoff renormalization group study Recently, Benyoussef and co- 

laborators have investigated the AAT using both mean-field treatment and Monte Carlo 
simulations ||21| , as well as finite-size scaling . In this paper we use a renormalization- 



group approach which, despite its simplicity, gives very accurate results for Ising-like 
models on the square lattice p3|,0, for the resistor network p4[ and directed percolation 



p5| . Motivated by these results we expect that, within its vahdity hmits, the present 
approach gives good estimates for the critical properties of the AAT, mainly for the phase 
boundaries. In the next section we define the model and obtain the recursion relations. 
In section III the results of our analysis are presented and compared with known results. 
Finally, in section IV we summarize our findings. 



II. THE MODEL AND THE RENORMALIZATION GROUP EQUATIONS 

The hamiltonian of the Ashkin- Teller model on a squate lattice is given by 

H = -Jl'^O-iaj - J2^TiTj - Ji^(TiTiajTj , (2.1) 
{ij} {ij) {ij) 

where cTj = ±1 and = ±1, Ji is the coupling between the spin variables a, J2 is the 
coupling between the spin variables r, and J4 represents the four-spin interaction that 
couples the two Ising systems. The sums runs over all first-neighbors pairs of sites. 

The model described by ( |2.1| ) will be studied by a real-space renormalization-group 
approach based on the self-dual graph shown in Fig. 1. It is well-known that the 



method is exact for classical spin models defined on the hierarchical lattice generated 
iterating the graph of Fig. 1, and produces a very good numerical approximation for the 
phase diagram of the square lattice. Since our approach does not allow any sublattice 
structure such as anti-ferromagnetic ordering we will require that the couplings satisfy 
the following ferromagnetic condition 

Ji + J2> 0, Ji + Ji> 0, J2 + > 0. (2.2) 

The renormalization is performed by the decimation of the spin variables associated 
with the sites 3 and 4 in Fig. 1. To simplify the calculations is convenient to introduce 
the transmissivity vector t = {ti,t2,t3) |26,2^ through 



^ ~ 1 + exp{-2Ki - 2K2) + exp{-2Ki - 2Ki) + exp(-2K2 - 2Ki) ' ^ ' ' 

^ ^ 1 + exp(-2iri - 2K2) - exp{-2Ki - 2K^) - exp(-2K2 - 2K4) 
^ ~ 1 + exp(-2fs:i - 27^2) + exp{-2Ki - 2K^ + exp(-2ir2 - 2/^4) ' 
_ 1 - exp(-2iri - 27^2) - exp(-2iri - 2K^) + exp(-2ir2 - 2K^) 
^ ~ 1 + exp(-2iri - 2K2) + exp(-2Js:i - 2K^) + exp(-2K2 - 2/^4) ' ^ ' ^ 
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where = Ji/ksT = pJi {i = 1, 2, 4). 

It can be shown that in the renormahzation of a series array of two bonds the vector 
t is simply given by the product of the corresponding vectors of the bonds ||26|j27[| . In a 



similar way, in the renormahzation of a parallel array of two bonds the same result holds 
for the dual vector t"^, whose components are given by 

= -—^ 2.6 

tf = ^-t^-t^ + t\ (2.8) 

The renormalization-group equations are obtained by requiring that the partition func- 
tion remains invariant after the decimation of the intermediate spins siting on the vertices 
3 and 4 represented in Fig. 1(a). The is achieved by writing 

eM-PK2 + K)= E E eM-Pni,2,3A) (2.9) 

{0-3, Ts} {u4,Ti} 

where T-C[ 2 and Tii, 2,3,4 are the cluster hamiltonians associated with Figs. 1(a) and 1(b), 
respectively, and K'f^ is an additive constant. 

It follows that the recursion relations for the renormahzation group thus defined are 
given by 

t[ = R(t\ + + titltl + tltl + 2tit2tl + 2titlh), (2.10) 
^2 = R{tl + + Ahtl + tjtl + 2tit2tl + 2t\t2h), (2.11) 
4 = R{tl + tl + t\flh + t\fl + 2titlh + 2t\t2h) (2. 12) 

where 

R = 2{l + t\ + t\ + tl + 2tl + 2tl + 2tl + 2titltl + 2tlt2tl + 2tlt%)-^ . (2.13) 

These equations completely determine the phase diagram for the AAT model, as well 
as the thermal critical exponents. Before we proceed with the analysis of the phase 
diagram, let us consider another set of variables defined trough 
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X = ^L±^, (2.14) 

r = tf, (2.15) 

Z = tf-tf. (2.16) 



Using the above introduced variables, the recursion relations given by Eqs. (|2.1CI| - |2TT^) 
become 

N^ No Nt 

X' = ^, r' = ^, Z'='^, (2.17) 

where 

A^i = 16X^X2 + 32X^Y + 16X3 + 48X2^^ + 16X2 - AXY'^Z^ - SXYZ"^ 

+12X -4Y^Z^ + 4:Z^, (2.18) 
N2 = IQX^Y + 16X^ + 64X3F - SX^YZ^ - 8X^Z^ - 16X1^^2 + 16^^ 

+16Y^ + YZ^ + Z^, (2.19) 
X3 = 48X2^ - 16X^^2^ - 32X^YZ - 32X^2^ + 32XZ + AY^Z^ 

+8YZ^ + 4:Z^, (2.20) 
D = IGX^Y + 16X^ + 32X3^=^ + 32X3 _ gx^YZ"^ + 24X^^2 - SXY'^Z'^ 

+24X^2 + 8Y* + IQY^ + YZ^ + + 8. (2.21) 

The new variable Z gives a measure of the anisotropy (essencially the difference be- 
tween K\ and -ft'2), and vanishes in the isotropic limit. We also note the Eqs. (|2.17| - |2!2l|) 
produce a symmetric flux under a reflection with respect to the Z = plane, in such a 
way that the regions Z > {K2 > Ki) and Z < {K2 < Ki) are isomorphous one to 
the other. For this reason we will restrict ourselves to discuss the subspace defined by 

0<X<1, 0<F<1, Z>0. (2.22) 

In the next section we will analyse the fixed points of our recursion relations, their 
respective thermal exponents and present the resulting phase diagram. 
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III. THE PHASE DIAGRAM 



The renormalization group expressed by Eqs. ( p.l0| )-( p.l2| ) reveals the existence of five 
phases separated by two-dimensional critical surfaces: (i) Paramagnetic P {< a >=< 
T >=< ar >= 0); (ii) Ferromagnetic F {< a >^ 0, < r >^ 0, < or >^ 0); (iii) 
Intermediate J (< cr >=< r >= 0, < ar >^ 0); (iv) Ferro-sigma F^- (< a >^ 0, < 
r >=< ar >= 0); and (v) Ferro-tau Fr (< r 0,< a >=< ctt >= 0). The resulting 
phase diagram in the space variables tf,t^ and is shown in Fig. 2, and reproduces 



qualitatively the results presented by |T^, |^ and |2^. The critical surfaces represent 
the boundaries of the domains of attraction of the trivial fixed points. On these surfaces 
we find nine (non-trivial) critical fixed points, which are presented in Table | along with 
their respective critical exponents uj- and determined from the scaling factor b and 
the relevant eigenvalues [^. The location of the non-trivial fixed points recovers all the 
available exact results for the square lattice. The fixed point Potts corresponds to the 
four-state Potts model and is completely unstable, with three relevants eigenvalues and 
two equal crossover exponents. All other non-trivial fixed points (Ji through Jg) are on the 
Ising universality class. Among these points, /i, I4 and I1 have two relevant eigenvalues. 
In particular, it is well-known that Ji is not in the Ising universality class 0; we expect 
that the same is true for I4 and Ij. The remaining non-trivial fixed points have a single 
relevant eigenvalue, and so they belong to the Ising universality class. As in the lAT, all 
phase transitions are second- or higher-order ones. 
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TABLE I. Non-trivial critical fixed points and critical exponents {vt,4>)- For the square 
lattice the exact values are fT = 2/3 and 1/7^ = 1 for the Potts and Ising universality class, 
respectively. 
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Fig. 3 shows the phase diagram in the variables X, Y and Z. As mentioned in the 
preceeding section, it suffices to consider the region Z > 0. This symmetry corresponds to 
a permutation of the a and r spin variables (Fo- and Fr phases) . In this diagram there is 
a self-dual plane which is invariant under the renormalization group transformation and 
is given by 

2X + Y + Z = 1. (3.1) 

This self-dual plane contains two lines of fixed points. One of these lines is the in- 
tersection of the self-dual and the Z = planes, corresponding to the lAT model and 
contains the F — P boundary. The other line is determined by the intersection of this 
self-dual plane with the other one given by 

2X + 4Y-Z = 2, Z>0, (3.2) 

and corresponds to the merging of the critical surfaces P — 1,1 — F^F — F„ and F„ — P. 

We believe that these results represent a very good approximation for the phase dia- 
gram of the AAT on the square lattice, since all the exact known results concerning the 
location of the critical frontiers were recovered. Unfortunately, as any Migdal-Kadanoff- 
like renormalization-group scheme, the present one fails to detect some special features 
of the universality classes, in particular the well-known result that the F — P boundary 
for the lAT is a line of varying critical exponent. In the following section we present our 
conclusions. 



IV. CONCLUSIONS 

We have used a real-space renormalization group approach to obtain the global phase 
diagram for the anisotropic ferromagnetic Ashkin-Teler model. Our findings can be seen 
either as an approximation for the AAT on the square lattice, or as an exact result for the 



hierarquical lattice generated by the scheme presented in Fig. 1 [p8| , p9| . In the former case, 
the resulting phase diagrams reproduces qualitatively well the structure of those obtained 
previously by other methods ||19|^I1. Since many exact known results for the square 
lattice are reproduced, we expect that the phase boundaries (critical surfaces) represent 
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a fairly good approximation for the exact ones. It also important to note that this simple 
approach gives non-trivial critical exponents, although it presents some discrepancies with 
respect to the universality class. In particular, the line of varying critical exponent is not 
reproduced. We believe that the results for the universality class can be improved if one 
considers sistematically larger clusters {b oo). However, such an enterprise is beyond 
the purpose of the present work, which was motivated by a relatively simple method to 
obtain reliable phase diagrams. 
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Figure Captions 



Fig. 1: Self-dual graph used for the generation of the hierarchical lattice considered in 
this work. The corresponding scaling factor is 6 = 2. 

Fig. 2: Phase diagram in space. 

Fig. 3: Phase diagram in [X, Y, Z) space. The Ft boundaries are obtained by a reflection 
of the Fg. boundaries on the Z — plane. 



12 



